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ABSTRACT 

We use a 3D photoionization modeling tool 
to study the morpho- kinematic properties of the 
Planetary Nebula (PN) NGC 3132. We show 
that it is possible to reproduce the low resolution 
observations (spectra and images) with an ellip- 
soidal shell. However, high resolution observa- 
tions, as those showing a density variation along 
the nebula and the [O iii]A5007 velocity profiles, 
definitively rule out this description. We show 
that a bipolar Diabolo shape with a 40° rotation 
of the symmetry axis relative to the line of sight 
successfully reproduces the observed images, as 
well as the high resolution observations. 



1. Introduction 

The determination of the three-dimensional matter dis- 
tribution in planetary nebulae is essential for a knowledge 
of the nebula ejection mechanism. A precise understand- 
ing of the nebula geometry which produces the observed 
morphology can also provide constraints on stellar evolu- 
tion theories. The geometry and the density distribution of 
PNe are generally derived from line imaging, assuming that 
a surface brightness enhancement corresponds to a density 
enhancement. After a first descriptive study of the mor- 
phology of PNe (Curtis 1918), more detailed studies indi- 
cate two different interpretations for the observed morphol- 
ogy of these objects: a) the different observed morphologies 
are due to different projections of a same common basic 
structure (Minkowski & Osterbrock 1960, Khromov & Ko- 
houtek 1968); b) the different types in which PNe can be 
classified (bipolar, elliptical, etc.) correspond to different 
evolutionary stages (Balick 1987). In this second approach, 
the morphology can also be related to the nebular and stel- 
lar properties, particularly the mass of the progenitor star 
(e.g. Peimbert & Torres-Peimbert 1983, Calvet & Peim- 
bert 1983, Stanghcllini, Corradi & Schwarz 1993, Corradi 
& Schwarz 1995). This point of view is corroborated by the 
studies of Zhang & Kwok (1998). These authors, assum- 
ing a three-dimensional ellipsoidal shell with density vari- 
ations and a complete absorption of the ionizing photons, 
which are reconverted in H recombination photons, produce 
simulated radio and Ha images which are compared to the 
observations. The visual appearence of PNe were also com- 
pared with models obtained for ionization-bounded prolate 
shells by Masson (1990). Assuming an ellipsoidal shell with 
constant thickness and that the ionizing radiation is geo- 
metrically diluted, Masson (1990) compared theoretical and 
observed Ha images. However, a given emission line image 
results from the sum of line intensities along the line of 
sight. The exact location of the regions which produce each 
line depends on the gas density distribution, but also on 



the ionization distribution, which, in turn, depends on the 
characteristics of the ionizing star and on the nebular gas 
abundance. Consequently, different emission lines provides 
information on different regions of the nebula, and only a 
consistent photoionization model applied to a realistic gas 
structure can give an accurate density distribution on the 
nebula. 

In this paper we show that the geometry of a PN derived 
by fitting the observed image of a line emission and/or a line 
intensity ratio is misleading, since different geometries can 
reproduce the same observational data. Other constraints, 
coupled with a detailed 3D photoionization model, must be 
used for better defining the geometry, as illustrated by the 
planetary nebula NGC 3132. 

Based on emission line imaging, an ellipsoidal geometry 
for NGC 3132 is assumed by various authors. Such a geom- 
etry is then suggested in a first approach. Studies similar to 
that of Masson (1990) using a prolate shell could explain the 
observed morphology of NGC 3132. Notice that the ellip- 
soidal shell model of Zhang & Kwok (1998) also reproduces 
the general trends of the observed Ha image of NGC 3132. 
Using a geometry defined by two concentric ellipsoids and a 
3D photoionization code, Biissgen, Diesch & Grewing (1990, 
hereafter BDG90| ) fit the observed [O m] and [N n] line im- 
ages, as well as the emission line spectra taken through four 
slits. 

Spatio- kinematic models for NGC 3132 were proposed by 
Sahu & Desai (1986, hereafter 3D86| ) to explain the ob- 
served asymmetric double-peaked [O iii]A5007 line profiles 
in five positions. [5D86 assert that an ellipsoidal model with 
velocity and density asymmetries reproduces the observed 
expansion velocities and correctly explains the nature of the 
profile asymmetries. However, in their paper, only the fit 
to the line profile at the central position is shown. Line 
profiles, as a powerful tool for studying the nebula geome- 
try, are discussed by Morisset, Gruenwald & Viegas (2000, 
hereafter MGV2000). 

As it will be shown below, all the models assuming an el- 
lipsoidal emitting shell lead to the same problems: neither 



can reproduce the low central density observed by SD86 and 
Juguet et al. (1988) nor the observed asymmetric double- 



peaked profiles in regions far from the center of the neb- 
ula. After a description of the available observations of 
NGC 3132 (§||), we show that ellipsoidal shell models do 
not reproduce all the available observations (§j|). A "Dia- 
bolo" nebula successfully explaining all the main observa- 
tional trends is presented in The conclusions are out- 
lined in Sj5l 



2. Observational data for NGC 3132 

A consistent photoionization model giving a realistic de- 
scription of a PN has to account for all the observational 
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data coming from different types of observation. In the fol- 
lowing, a summary of the data available for NGC 3132 is 
presented. 

The logarithm of the observed H/3 flux ranges from -10.49 
flPerek 1971| ) to -10.20 ( |Pottasch et al. 1977|). The extinc- 
tion E(B-V) is of the order of 0.1 ( [Pottasch et al. 1977 , 
Mendez 1978, Feibelman 1982, Gathier, Pottasch & Pel 
1986). 

The calculated values for the Zanstra He n temperature 
of the ionizing star of NGC 3132 ranges from 73 000 K (de 
Freitas Pacheco , Codina & Viadana 1986) to 110 000 K 
(|Pottasch 1996; ). Distance determinations indicate values 
from ^0.51 kpc ( |de Freitas Pacheco et al. 1986| , Pottasch 
1996) to 1.63 kpc (Torres-Peimbert & Peimbert 1977). The 
ionizing star luminosity available in the literature is 72 Lq 
(Mendez 1978) or - 125 L ( |Pottasch 198^ ). 

HST/WFPC2 images for Ha, [O i]A6300, [O m]A5007, 
[N n]A6583, and [S n]AA6717,6731 are available at the 
CADC server from the Trauger's proposal # 622 ifj. Lower 
resolution images and contours were obtained by Juguet et 
al. (jjH| |BDG90|, [Pascoli (1990)| , pondal et al. (1981)| , 
and Schwarz et al. (1992). 

Spectroscopic observations of NGC 3132 can be found 
in Aller & Faulkner (1964), Kaler (1976), Torres-Peimbert 



a geometrical ellipsoidal shell. 



3. The ellipsoidal model 

The 3D photoionization code used in this paper is de- 
scribed in Gruenwald, Viegas & Broguiere (1997). The out- 
put of the code is treated by IDL tools performing rotations, 
projections, and velocity profiles. A full description of these 
tools is given by MGV2000. 

In order to fit the observed morphology of NGC 3132 
([O iii]A5007 and [N n]A6583 line image s), as we ll as the 
emission line spectra taken through 4 slits, BDG90 proposed 
a model obtained by a 3D photoionization code. The geom- 
etry and input parameters proposed by BDG9C are taken as 
a starting point for our 3D model of NGC 313 2. Elements 
heavier than Ne were not included in BDG90 model. For 
S, CI, and Ar we assumed a solar abundance ( prevesse fe 



Anders 1989), while Mg, S, CI, and Fe are depleted by one 



& Peir pbert (1977) 
14.7 km s -1 



and BDG90 



Expansion velocities of 
for [O m]A5007 and 21.0 km s" 1 for [O n]AA- 



hundrcd (Stasinska & Tylenda 1986). Thus, the gas chemi- 
cal abundance, in number, is: H: 1.0, He: 0.126, C: 7.1(-4), 
N: 2(-4), O: 6(-4), Ne: 8.2(-5), Mg: 3.8(-7), Si: 3.55(-7), S: 
1.62(-5), CI: 3.16(-7), Ar: 3.6(-6), Fe: 4.7(-7). The central 
ionizing radiation is a black-body of 90 000 K and 150 L Q . 
The model consists of two prolate ellipsoids defining a shell 
and an inner cavity (the outer zone being empty). The in- 
ner ellipsoid has a semi-minor axis of 1.27 10 17 cm and a 



3726,3729 were obtained by Meatheringham, Wood &Faulkner se mi- ma jor axis of 2.32 10 17 cm, while the outer one has 
(1988), for slits along the major axis. High resolution spec- 
tra for [O iii]A5007 in five positions were obtained by SD86 , 
who also deduced a velocity of 14 km s _1 at the central 
position. The line profiles are all asymmetrical and double 
peaked. Since the observed red peak is more intense to- 
wards the center, the asymmetry can not be due to local 
absorption. 

Electronic densities are generally obtained from the [O n] 
or [S n] line intensity ratios corresponding to intensities in- 
tegrated on a slit 



Torres-Peimbert & Peimbert (1977 



tained 1000 cm , while Meatheringham et al. (1988 



ob- 
and 

Stanghellini & Kaler (1989) give, respectively, 600 cm~ 3 
and 430 cm~ 3 . More detailed observations show a vari- 
ation of the electronic density measured in different posi- 
tions. Density determination by Juguet et al. (1988)| , from 
the [S n] AA6717,6731 line intensity ratio, shows a double- 
peaked distribution along the NS direction, reaching up to 
1300 cm" 3 at the periphery of the nebula and decreasing 
to 300 cm -3 at the central position. Recent observations 
have confirmed this density distribution trend (Monteiro, 
Gruenwald & de Souza, in preparation). It will be shown 
that the decrease in density at the central region is the key 
observation that rules out the description of NGC 3132 as 
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1.96 10 17 cm and 2.45 10 17 cm, respectively. The density in 
the inner cavity is constant (464 cm -3 ). In BDG90, the 
shell density varies along the three directions. Here only 
the most important density gradient (along the major axis 
Z) is assumed. It varies from ~ 1200 cm -3 at one pole to 
~ 900 cm -3 at the opposite pole. This gradient is needed 
in order to reproduce the observed asymmetrical brightness 
enhancement shown by the optical images. 

The calculations are performed assuming the on-the-spot 
approximation (see MGV2000). Initially the cube contain- 
ing 1/8 the nebula is divided in 36 3 cells. In order to im- 
prove the numerical resolution, some of these cells are sub- 
divided, resulting in 133 000 cells. Once the calculations 
converged, the whole nebula is recovered (see MGV2000). 
Following [BDG9C| , the images are obtained after a 7°. 5 rota- 
tion around the X axis (perpendicular to the line of sight). 



3.1. Matter- and Radiation-bound models 

The results obtained by BDG90] for the line intensities 
are not reproduced by our ellipsoidal model. First, the ge- 
ometrical thickness of the shell, adopted by these authors, 
is not lar ge enoug h to absorb all the photons. Our model 
with the BDG90| geometry leads to a matter-bound (M- 
bound) nebula. The model is then unable to reproduce the 
observed [N i] AA5198,5200 strong lines, characteristic of a 
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radiation bound nebula. Notice that the measured inten- 
sity for [O i]A6300, - 0.3 H/3 (Torres-Peimbert & Peim- 
bert 1977), also indicates that the nebula must be mainly 
radiation-bound (R-bound). Although the details of the 
BDG90I code are not available in their paper, the difference 



with our results may come from the optical depth calcu- 
lation, which is the key-parameter determining the model. 
Models with a low number of cells (which implies too large 
cells) tend to overestimate the optical depth. In this case, 
for the same geometrical thickness, the model would be R- 
bound ins tead of M -bound. Our model has 40 times more 
cells than BDG9C. This can explain why, using the same 



input parameters, our model is M-bound while theirs is R- 
bound. In order to have a R-bound shell, reproducing the 
observed low-ionization emission lines, the outer ellipsoid 
size must be increased by at least 20%. A second reason f or 
the discrepancies between our model and that of BDG9C is 
the lack of elements heavier than Ne in their code, which 
leads to an underestimation of the cooling processes. In or- 
der to test this effect, a heavy-element free-BDG90 model 
was built. The [O m]A5007 and [N n]A6583 line intensities 
(relative to H/3) increase, respectively, by 30% an d 16%. 
This can partly explain the differences with BDG9C results, 
since the optical depth effect, discussed above, must be also 
important. 

In the following, we will refer to the matter-bound model 
using the dimensions of BDG90 as the M-BDG90 model 
and to the 20% "extended" radiation-bound model as the 
R- BDG90| m odel. The results of both models are presented 
in Table |1] As expected, the stronger differences appear 
for the low ionization lines ([N i]AA5198 ,5200 and [O i]A- 
6300). Nevertheless, even the R- BDG90| mod el is far from 
reproducing the results obtained by BDG90, as discussed 
below. 



EDITOR: PLACE TABLE g HERE. 



3.2. Ellipsoidal model results 

The superposition of the [N n]A6583 HS T image and the 
theoretical isophotes obtained from the R- BDG90 model is 
shown in Fig. [j]. The shape is well reproduced, as well 
as the variation of the surface brightness along the major 
axis. To adjust the modeled shape and dimensions to the 
observed ones, a reduction of 4% in the distance used by 
BDG90 (670 pc) is necessary. Note that this new value for 
the distance is in the range of the calculated distances to 
NGC 3132 (§2) 

For a comparison of our results with the observations of 
Juguet et al. (1988), the theoretical N-S density distribu- 



model assuming an ellipsoidal geometry (dashed line) does 
not reproduce the observed central decrease in density. 

In order to obtain the theoretical [O m]A5007 emis- 
sion line profiles, a radial velocity law V = a.r/r + (3.r 
is used. Three possible sets of a and (3 values are as- 
sumed: a = 14.7,9.1,0 km s^ 1 and (3 = 0, 9.1/r , 23.8/r 
km s _1 cm -1 , respectively. The normalization factor, r a 
= 3 10 17 cm, is of the order of the size of the nebula. The 
< V >5oo7 and < V >3726 emission line mean velocities for 
the whole nebula (as defined in MGV2000) are given in Ta- 
ble ||. For the three sets of values defining the velocity law, 
the parameters were chosen in order to give < V >5oo7 = 
14.7 km s -1 , as obtained from the observed profile reported 
by Meathcringham ct al. (1988)| . The corresponding theo- 



tion is calculated from the theoretical [S n]AA6717,6731 ra- 
tio (using the Alexander & Balick 1997 fit) for a 4" width 
slit, centered on the nebula. As seen from Fig. 0, this 



retical < V >3726 mean velocity depends on the parameters 
a and (3 and is always slightly lower than the observed value 
(21 km s _1 ). The five positions at which the [ O iii]A5007 
line profiles were observed are shown in [5D86 (their Fig. 
la) superposed on an U-band image (Kohoutek & Lausten 
1977). The calculated velocity profiles for [O iii]A5007 are 
shown in Fig. |^. Since the three sets of parameters of 
the adopted velocity law give similar line profiles, only the 
results for the a = 9.1 km s _1 and (3 = 9.1/r G km s _1 
cm -1 law is shown. In the upper- right panel the aperture is 
shifted from the center by 2.8" NE. This shift corresponds 
to the pointing accuracy of the telescope used by 5D86| and 
it allows to check if the effect of the pointing precision on 
the profile is important. The central double peak is re- 
produced. However, the observed asymmetry can only be 
obtained if the nebula is tilted by about 30°, whereas the 
results shown in Fig. ^ were obtained assuming a rotation 
angle of 7°. 5, as proposed by BDG90. The observed dou- 
ble peaks at the four external positions are not reproduced. 
Note that the observed profiles were taken by SD86 at posi- 
tions where the U band emission is strongest. In the U band, 
the nebular emission is dominated by the [O n]A3726 line, 
so the observations were taken at the extreme outer part 
of the nebula. Thus, the radial velocity at these positions 
is quasi-perpendicular to the line of sight and an apparent 
velocity of ~ 10 km s _1 would require an unrealistic high 
radial velocity. 

EDITOR: PLACE TABLE | HERE. 

Finally, a simple way to verify the inadequacy of the ellip- 
soidal model in reproducing the density gradient and the low 
ionization lines is the following: let us consider a dense shell 
(~ 1000 cm~ 3 ) with a geometrical thickness t and radius R, 
around a less dense cavity (~ 300 cm~ 3 ). The line emissivity 
is proportional to the square of the density. Thus, in order 
to reproduce the observed density distribution, the main 
contribution to the central emission must come from the 
low density gas, implying t < 0.02R. Such a thin shell will 
have a total emissivity far below the observed one, and will 
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be matter-bound, not producing the low ionization lines. A 
radiation-bound model, necessary to explain the low ioniza- 
tion lines cannot reproduce the density gradient as shown 
above. This rules out the description of NGC 3132 as a 
closed ellipsoidal shell. An ellipsoidal shell with a hole just 
aligned to the central line of sight could be proposed, but 
such an ad hoc geometry seems unrealistic. 

From all the results presented above we conclude that an 
ellipsoidal model for NGC 3132 is not appropriate to explain 
the observational data. Thus, a new model is proposed be- 
low. 



4. The Diabolo model 

PNe imaging suggests that several objects can have a "Di- 
abolo" shape (e.g. Hour-Glass nebula and NGC 2346). In 
order to reproduce this shape, we adopt a two-zone gas dis- 
tribution, illustrated in Fig. ^: a) a dense zone (with a 
density of 1300 cm -3 ) defined by the intersection of two 
spheres of the same radius (3. 10 17 cm) but with centers 
shifted by 10 17 cm, defining the "Diabolo-shape" , and b) 
a low-density gas, 300 cm~ 3 , filling the rest of the nebula. 
The shadowing in the figure is just used to give a better 
3D visualization and does not correspond to a density vari- 
ation. The model was calculated for 1/8 of the nebula and 
the whole nebula was reconstructed in a 10 6 cells cube. We 
keep the chemical abundances and ionizing star character- 
istics of BDG9C model (§||). No attempt to adjust these 
parameters were made in order to improve the fit of the 
emission line ratios. In fact, the scope of this paper is just 
to show that a strong constraint on the geometry is the 
density distribution. Once this distribution is reproduced, 
other observed features, unexplained by previous models, 
can also be explained. 



4.1. The Diabolo model results 

The effect of varying the angle of view of the nebula is 
shown in Fig. ||. Practically all the observed regular mor- 
phologies of PNe are reproduced: from a clearly butterfly 
image (axis of symmetry perpendicular to the line of sight) 
to a well round nebula (axis of symmetry aligned to the 
line of sight) through an ellipsoidal one. With an appro- 
priate orientation of 40° relative to the line of sight, the 
observed optical shape of NGC 3132 (ellipsoidal shape) can 
be reproduced. As said above (§1), such feature can also 
be reproduced by the ellipsoidal shell model of Zhang & 
Kwok (1998), as well as by a model similar to that of Mas- 
son (1990). However, these models are not self-consistent 
with respect to the ionization balance and, as will be shown 
below, other constraints will rule out an ellipsoidal model 
for this nebula. 

With a Diabolo morphology the contribution from the 



dense region to the central emission is reduced, helping to 
solve the close ellipsoidal shell problems. At the same time, 
the projected velocity of the gas in the external parts of 
the nebula will be high enough, reproducing the observed 
double peak, as shown below. 

Regarding the emission line ratios, the results of the Dia- 
bolo model (fourth column of Table [l]) are similar to those 
obtained from the R-BDG90 model. In addition, HST line 



images in H/3, [O m]A5007, [N n]A6583 and [O i]A6300 are 
shown Fig. ^ (four upper panels) and can be compared to 
the corresponding Diabolo model images (four lower pan- 
els). The observed ellipsoidal shape is well reproduced, as 
well as the ionization stratification. The brightness asym- 
metry along the major axis is not reproduced by the adopted 
axi-symmetrical model of the nebula. Such asymmetry 
could be explained either by a density gradient or an ion- 
ization source that is not in the geometrical center of the 
Diabolo. 

The electronic density distribution obtained from the 
[S n]AA6717,6731 ratio is also shown in Fig. || (solid line). 
The observed decrease of the density (§2) in the central part 
of the nebula is reproduced. This and the line profiles dis- 
cussed below are the main improvements obtained by the 
Diabolo shape. 

The same velocity laws applied for the R-BDG90 model 
are used (§3.2) with r D 



5 10 17 cm. The three different 



[O m]A5007 profiles obtained with the three sets of param- 
eters of the velocity law are shown in Fig. |^ and can be 
compared to the observations (3D86). The profiles show 
the double peaks at the six positions, a feature not repro- 
duced by the ellipsoidal shell model. The asymmetry is well 
reproduced at the four external positions. The asymmetry 
at the central position is not reproduced, as long as an axi- 
symmetric nebula is used, even considering a shift of 2.8" as 
shown in the upper right panel. Notice that at the central 
position only the low density gas contributes to the emis- 
sion. Thus, any asymmetry in the velocity field and/or in 
the low density distribution will lead to the observed asym- 
metry, requiring no change in the higher density distribution 
associated to the Diabolo shape. 

The three sets of parameters for the velocity law were cho- 
sen in order to obtain the same < V >5oo7 (14.7 km s" 1 , 
Table ||; the results for < V >3726 remain the same as for 
the R-BDG90 models); however, for a — (dashed lines), 
the peaks of the observed profiles at the central part of the 
nebula occur for a lower velocity (Fig. |?], upper panels). 
Furthermore, in this case, the velocity law does not repro- 
duce the observed profiles since one of the peaks is almost 
inexistent. The results obtained with the other two sets 
of parameters (solid and dot-dashed lines) are very similar. 
Nevertheless, for a / and [5 ^ 0, < V >3726 is higher 
than < V > 5007 5 approaching the observed result (Table |^) 



of Meatheringham et al. (1988) 



5 



5. Concluding remarks 

We have shown in this paper that the lack of high reso- 
lution observations (spatial and spectroscopic) may lead to 
a wrong conclusion about the morphology of a PN. 

A new modeling tool based on a 3D photoionization code 
is used to consistently reproduce all the available observa- 
tions of NGC 3132. Using an ellipsoidal geometry, all the 
low resolution observations (flux, emission line spectrum, 
imaging) could be reproduced. However, neither the den- 
sity variations indicated by the [S n]AA6717,6731 ratio, nor 
the [O iii]A5007 velocity profiles could be explained by such 
a model. To account for these observations, a drastic change 
in the geometry of the nebula is necessary. The proposed 
Diabolo shape offers the solution to the density distribu- 
tion, as well as to the asymmetric double-peaked emission 
line profiles, while also explaining the low resolution obser- 
vations. 

A Diabolo shape may reproduce the observed bipolar 
morphology of many PNe. In fact, Bryce, Balick & Meaburn 
(1994) suggested that NGC 6720 is not a spherical expand- 
ing shell as often assumed, but must be bipolar in nature, 
since the observed [Nil] emission line profiles are split in two 
components. A bipolar shape emerges from studies of the 
formation and evolution of PNe, in particular those related 
to the interacting wind model (Kwok, Purton, & Fitzger- 
ald 1978). Recent hydrodynamical models (e.g. Dwarkadas, 
Chevalier & Blondin 1996; [Mellema 1997] ) suggest that a Di- 
abolo type shape can be produced in the accepted scenario 
of PNe formation (post AGB envelope ejection with wind 
interaction). In fact, recent high resolution images taken 
with the Hubble Space Telescope show that PNe with this 
type of structure is not uncommon. Some nice examples are 
the Hour-Glass nebula, NGC 2346, NGC 6537 and Hb-5. 

No attempt was made to fine tuning the emission line 
ratios. This can be achieved through changes in the ioniz- 
ing star characteristics, in the chemical abundances, as well 
as through small adjustments in the density distribution, 
which will not change the main conclusions of this paper, in 
particular the Diabolo shape. Furthermore, in the models 
presented here, the ionizing source is assumed to be at the 
center of symmetry of the nebula. However, there are evi- 
dences that the star should be shifted by ~ 1.7" from the 
geo metrica l center of the nebula (Kohoutek & Lausten 1977 
and 5D86 ). A model accounting for this out-of-center ion- 
izing star may still improve the results, leading to a better 
agreement with the observations regarding the brightness 
asymmetry and the emission line profile at the central po- 
sition. 

We have shown here the importance of a 3D photoion- 
ization model, combined with imaging results and also with 
high resolution observations, to impose relevant constraints 
on the PN geometry. The misclassification of NGC 3132 as 
an ellipsoidal PN instead of a bipolar one casts doubts on 



the statistical results correlating the morphology and other 
PNe properties, such as abundances or progenitor type. 
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Fig. 1. — Superposition of the HST and the theoretical im- 
ages of [N n]A6583. The theoretical results correspond to 
the R-|BDG90| model. 



Fig. 2.— N<->S variation of [S n] AA6717,6731 density from 



the R- BDG90 model (dashed line) and from the Diabolo 
model (solid line). 

Fig. 3. — [O iii]A5007 velocity profiles obtained from the 
R- BDG90| model. The positions and the size aperture are 
the same as in |3D86 (their Fig. 2). In the upper-right panel 
the aperture is shifted from the center by 2.8" (see text). 

Fig. 4. — Gas distribution for the Diabolo model. Only the 
denser zone (1300 cm" 3 ) is shown. 

Fig. 5. — H/3 images from the Diabolo model. From the 
upper-left to the lower-right panel the angle between the 
axis of symmetry and the sky plane increases by 10° from 
0° to 80°. 

Fig. 6. — The top panels correspond to HST images of 
NGC 3132 while the botton panels to the Diabolo model. 

Fig. 7. — [O iii]A5007 velocity profiles for the Diabolo model 
using the velocity law V = a.f/r + f3.r, with ro = 310 17 cm, 
and (a, 0) = (14.7, 0), (9.1, 9.1/r ), and (0, 23.8/r G ) are 
shown, respectively, by dashed, solid, and dot-dashed lines. 



This piano tables was prepared with the AAS iATf^X macros v4.0. 
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Table 1: Emission line intensities relative to H/3 



Line A. 


M-lRnciqol 

IVlTlJ 1 J V.I C/ Ul 




UlCl UU1U 




fl 099 


u.uio 


U.UIO 




U.UOtJ 




U.UOJ 


111J AOUU i 


O.U40 


/i 7nn 
4. / uu 


o.oyu 


[N t A5200 


0.010 


0.159 


0.200 


He IA5876 


0.175 


0.189 


0.190 


[0 i]A6300 


0.037 


0.320 


0.353 


[N n]A6583 


2.792 


5.090 


5.337 


[S n]A6717 


0.193 


0.689 


0.825 


[S n]A6731 


0.236 


0.745 


0.846 


Ha 


2.904 


2.915 


2.917 


H/3 1 


6.553 


9.651 


9.050 


In 10~ 11 erg cm 


2 s 1 , for a distance of 670 pc. 




t^: p. i 

Table 2: Mean line velocities for the adopted velocity law °' ' 


a a 0" 


< V >5007 Q 


< V >3726 a 




14.7 0. 


14.7 


14.7 




9.1 9.1 


14.7 


15.9 


1000 


0. 23.8 


14.7 


17.8 




a _i 
m km s 






800 



in km s _1 /r , with r D = 310 17 cm for the R-BDG90 model and r = 510 17 cm 
for the Diabolo model 
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